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NASTRAN MODEL OF A LARGE FLEXIBLE SWING-WING BOMBER 


Volume V; NASTRAN Model Development— Fairing Structure 

W. D. Mock and R. A. Latham 
Rockwell International Corporation 
Los Angeles, California 


SUNMARY 


This report describes the development and validation of tfie NASTRAN 
model of the B-1 aircraft 2 (A/C-2) fairin? structure. The development of 
this model completes the seven substructures defined in the Airloads Research 
Study NASTRAN Model Plans. Subsequently, these seven NASTRAN substructure 
models will be joined to form the NASTRAN model representation of the total 
aircraft structure. The intent is to utilize the asse.mbled MASTRAN model 
computed stif&iess matrix in conjunction with the FLEXST.AB program for 
aeroelastic analysis. The application of these advanced programs on a large, 
flexible aircraft that has accumulated significant flight -test data will add 
to the technology base for future transport aircraft. 

During this contract phase, the NASTRAN model plan for the fairing 
structure was e.xpanded in detail to generate the NASTRAN model of this sub- 
structure. The grid point coordinates, element definitions, material 
properties, and sizing data for each element were specified. 

The fairing model was throughly checked out for continuity, connectivity, 
and constraints. The substructure was processed for structural influence 
coefficients (SIC) point loadings to determine the deflection characteristics 
of the fairing model. Finally, a demonstration and validation processing of 
this substructure was accoirplished using the .NASTRAN finite -element program 
installed at the NASA DFRC facility. The bulk data deck, stiffness matrices, 
and SIC output data were delivered to .NASA DFRC. 


LNTRODUCTION 


A/C -2 (figure 1) is employed in the .Airloads Survey Flight Test program. 
This aircraft has undergone extensive gro'ond testing to calibrate the strain 




gages utilized in the airloads surv'ey. The aircraft provides a reasonable 
sijitulation of a future transport aircraft since it enploys the large, flexible 
structure (figure 2) envisioned in future transport designs. 

The airloads data gathered during the flight -test program can be utilized 
in the evaluation of NASA computer programs recently developed to enhance 
the analytical techniques of predicting aeroelastic response of large, 
flexible aircraft. These analytical techniques include computerized structural 
analysis programs such as NASTRA.N and FLEXSTAB. 

Since the B-1 development program involves all e.xperimental tests needed 
to correlate the analytical pi'edictions with actual measured results, detailed 
plans for constructing a XASTRAN structural model of the B-1 airframe suitable 
for use on the NASA/DFRC Cyber computer were initiated. This model is of 
minimum conqjlexity to give satisfactory flexibility characteristics for the 
FLEXSTAB aeroelastic analysis. Included in this model are the control 
surfaces, control system stiffiiess, and secondar>’ leading edge and trailing 
edge structure. During this contract phase, detailed plans for constructing 
a MASTRAN model of the fairing substructure were implemented. Grid point 
coordinates for this substructure were coded for each element, and the material 
properties and sizing data were specified. The bulk data were throughly 
checked using interactive graphics techniques. The data were evaluated for 
continuity, connectivity, and constraints. In addition, the SIC point load- 
ings were applied to compute the deflections at selected locations. A 
demonstration and validation processing of the N.ASTR%V model substructure was 
accomplished using the XASTRAN finite-element program installed or. the .NASA 
DFRC Cyber computer. 


AIRCRAFT DESCRIPTION 


The B-1 aircraft is a prototype long-range supersonic bomber with the 
capability of high-speed flight at low altitude. Configuration dimensions 
and general arrangement are presented in figure A-1. The aircraft utilizes 
a blended ;ving-body concept with variable -sweep wings, a single vertical 
stabilizer with a three-section (upper, intermediate, and lower) rudder, 
and horizontal stabilizers which operate independently to provide both pitch 
and roll control. The variable -sweep (15 to 67.5 degrees) wing, equipped with 
slats, spoilers (which also function as speed brakes), and flaps, provides the 
aircraft with a highly versatile operating envelope. Canted vanes, mounted 
on each side of the forward fuselage, are part of the structural mode control 
system which reduces structural bending oscillations in the vertical and 
lateral axes. 
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The aircraft is powered by rour \T10l-GE-100 dual -rotor augnented 
turbofan engines in the 30, 000 -pound- thrust class. The engines are mounted 
in twin nacelles below the wing, approximately at the left and right wing 
pivot points. For supersonic speeds, an air induction control system varies 
the internal geometry of the nacelle inlet ducts to maintain the required 
airflow to the engines for all flight conditions. 


FAIRING 


The total fairing substructure consists of the ov'^enving fairings sho\vn 
ui figure 3 and the undenving fairings shown in figure 3. The fairing 
support structure is shown in figure 4. 

The overwing fairing is comprised of the upper pivot fairings, the 
forward intermediate fairing, and the ovenving movable fairing. The upper 
pivot fairings are above the wing pivot fitting and cover t.he region from 
fuselage stations 875 to 1036 and butt lines 119 to 188. These fairings are 
laminated fiberglass panels which are preloaded against the wing surface. 

Tlie forward intermediate fairing, aft of the pivot fairings, is a sandwich 
panel with an aluminum core and fiberglass cover. This panel is cantilevered 
from the fuselage. The overwing movable fairing is a full -depth fiberglass 
honeycomb panel which extends aft to approximately fuselage station 1140. This 
movable fairing has a hinge support mounted on the fo^^^^ard intermediate fairing 
panel. The movable fairing is actuated by a track/ trolley arrangenent mounted 
on the movable fairing with support connection to the wing inboard trailing 
surface structure. 

riie underwing fairings are comprised of the lower pivot fairing and the 
intermediate and aft underwing fairing panels. The lower pivot fairing is 
below the wing pivot fitting and covers the region between fuselage stations 
875 and 993 from butt lines 119 to 188. This fairing segment is constructed 
of laminated fiberglass and is preloaded against the wing stixicture. The 
intermediate and aft panels extend over the top of the nacelle structure and 
are supported by a series of linkages mounted on the nacelle structure. These 
linkages enable the panels to be actuated up or down during the wingsweep 
operations. The lower intermediate fairing panel is a machined aluminum plate. 
The lower aft panel is a sandwich panel with aluminum core and fiberglass face 
sheets . 

Figures 6 through 8 are photographs showing the overwing fairing during 
the wing sweep operation. The pivot attach region of the support structure is 
shown in figure 9, viewed from the right-hand side of the aircraft, looking 









Figure 6 


Fairing position at forward sx^eep 



o 



Figure 7. - Fairing position at intermediate sweep 






Figure 9. - Fairing support structure to wing attachment. 




foTOard. Figure 10 shows a partial view of the under.ving fairLug linkage 
mechanism. 


NASTR.AN MODELS 


The detailed plans for the finite -element modeling of the A/C -2 
structure intended for use with the NASA COSMIC release of NASTRAN level 16.0 
on the NASA DFRC Cyber computer constrains th.e model to the minimum complexity 
to give satisfactory flexibility characteristics for FLEXST.AB aeroelastic 
analysis. 

The NASTRAN model plans specify seven substructures consisting of the 
following: 

(1) Horizontal stabilizer; leading edge, and trailing edge 

(2) Vertical stabilizer; leading edge, and rudders 

(3) Nacelle structure 

(4) Wing outer panel, flaps, slats, and outboard transition ribs 

(5) Forvv'ard fuselage structure 

(6) Aft fuselage structure, wing carr^/- through structure (WCTS) , and 
inboard transition lugs 

(7) Overwing and underlining fairings 

In addition to modeling the A/C -2 airframe structure to represent the 
flexibility characteristics, the model was designed to provide stress data 
at the airload survey strain gage locations for each component, in these 
regions, the model complexity was increased to provide the desired accuracy. 

In some regions, the complexity was dictated by the NASTRAN aspect -ratio 
constraints. During this contract phase, the NASTRAN model plans for the 
fairing structure were implemented to generate the NASTRAN model for this 
substructure. The description of this model, which was demonstrated and 
validated on the NASA DFRC Cyber computer system, is presented herein. 
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Fairing NASm^A' Model 


The NASTRAN fairing model is configured to be representative of the 
A/C-2 structure. The aluminum structure which forms the leading edges of the 
upper and lower pivot fairings beti>reen fuselage stations 863 and 875.5 has beer- 
idealized using the NASTRAN plate elements CQUAD2 and CTRIA2. The fuselage 
station 875.5 aluminum bulkhead has been idealized as a series of conrods and 
shear pjatels. The laminate fiberglass sections and aluminum region of the 
fairing structure have been idealized by the NASTRAN plate elements CQUAD2 
:md CTRIA2, wich include bending, n«r±»rane, and transverse shear effects. The 
sandwich contruction regions have been modeled utilizing the NASTR.AN CQU.ADl 
:and CTRIAl sanchvich elements, which include bending, membrane, and transverse 
shear effect. The support structures for the pivot fairings are represented by 
C3<\R elisments with titanium material properties , except for the upper beam at 
fuselage station 992 and the lower outpoard beam extending n-om fuselage 
stations 884 to 944, which are made of steel. The remaining fairing support 
structure is represented by GEAR elements using aluminum material properties. 
The coniiections of the various fairing segments at the slip-slide joints are 
represented by the .NASTRAN CELASl elements. 

.A summary of the number of grid points and element t^'pes utilized for 
the model idealization of the fairing structure is presented in table I. 

The model diagrams of the fairi’^g substructure which define the eleiwn.t 
types, element ID numbers, and grid numbering systems are presented in figures 
C-1 through C-8. 

The fairing is constrained vertically along the leading edge at fuselage 
station 863 and along the fuselage interface from fuselage stations 865 
through 1096. The fairing is also constrained in the inboard -outboard direc- 
tion along the fuselage interface. The beam at fuselage station 944 is con- 
strained vertically at butt line 148.0. The overwing movable fairing is 
•restrained vertically at fuselage station 1108 and at butt line 154 to 
■represent the fairing -to -wing tie. The unden>fing fonvrard fairing panel is con- 
strained vertically at the four linkage points at fuselage stations 1142 and 
1139, ITie underwing aft panel is constrained in three directions at fuselage 
station 1243.5. 

The -Airloads Research Study NASTRAiN model was thoroughly chec.ked out for 
continuity, connectivity, and constraints using interactive graphics tech- 
niques. This model was then processed for the loading applied at each SIC 
point (table C-1 and figures C-9 and C-10) , with fairing structure appropri- 
ately supported. The deflections computed for these SIC loadings are shoun 
in figures C-11 through C-15. 
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TABLE I. - ARS NASTRAN MODEL STATISTICS 


Description 

of 

substructure 

.NASTRAN nradel elements 

No. 

of grids 

Rods 

Bars 


Sandwich 

plates 

Plates 

Scalar 

Springs 

Fairings 

504 

10 

167 

3 

64 

120 

“6 


Element 

N.ASTRAiN nomenclature 

Rod 

= COK'ROD 

Bar 

= CBAR 

Shear panel 

= CSHEAR 

Sandwich 

* CQUADl and CTRIAl 

Plate 

= CQU.AD2 and CTRIA2 

Scalar spring 

= CEL.AS1 


BULK DATA 


The NASTRAN model coordinates, sicing, material properties, and loading 
data are presented in the NASTRAN program input format . Since these data are 
identified by column numbers excerpts from the NASTRAN User's Manual which 
define the contents of the card columns for each card tvpe are included. The 
format of the sorted bulk data for each element t>T5e is presented on pages 25 
through 45 . This format is applicable to the NASTR.AN model bulk data pre- 
sented on pages 72 through 93. 
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FIGURES USING ENGINEERING UNITS 
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Figure A-1. - General arranganent - RDT^E A/C-1 and-2 
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Appendix B 

NASTRAN NDDEL BULK DATA FORMAT 




BULK DATA DECK 


Input Data Card CBAR Simple Beam Element Connection 

Oescn'ption : Defines a simple beam element (BAR) of the structural model. 

Format and Example : 


1 2 3 A 5 6 7 8 9 10 


nMii 


PIO 

GA 

GB 

XI ,G0 

X2 

X3 

F 

abc 

CBAR 

2 

39 

7 

3 

13 



2 

123 

+bc 

PA 

PB 

Z1A 

Z2A 

Z3A 

ZIB 

Z2B 

Z3B 

WBB 

♦23 


513- 








— 1 


Field 

EID 

PIO 

GA.G8 

I X1.X2.X3 


GO 


Contents 

Unique element identification number (Integer > 0) 

Identification number of a PBAR property carp (Default is EID unless 3AR0R card 
has nonzero entry in field 3) (Integer > 0 or blank*) 

Grid point identification numbers of connection points (Integer > 0; GA * GB) 

Components of vector v. at end a, (figure 1(a) on page 1.3-15) measured at end a, 
parallel to the component: of the displacement coordinate system for GA, to deter- 
mine (with the vector from end a to end b) the orientation of the element coordinate 
system for the bar element (Real, X1^ + X2* + X3^ > 0 or blank*, see below). 

Grid point identification number to optionally supply XI, X2, X3 (integer > 0 or 
blank*) (see below) 


F 


I PA,PB 


Flag to specify the nature of fields 6-8 as follows: 
6 7 8 


F = blank* 




F * 1 

XI 

X2 

X3 

F . 2 

GO 

blank/0 

b1ank/0 


Pin flags for bar ends a and b , respectively, that are used to insure that the bar 
cannot resist a force or moment corresponding to the pin flag at that respective 
end of the bar. (Up to 5 of the unique digits 1-6 anywhere in the field with no 
imbedded blanks; integer > 0) (These degree of freedom qodes refer to the element 
forces and not global forces. The bar must have stiffness associated with the 
pin flag. For example, if pin flag < is specified, the bar must have a value for 
J, the torsional constant.) 


Z1A,Z2A,Z3A Components of offset vectors Wj and W(j, respectively, (see figure 1(a), page 
Z18,Z2B,Z38 1.3-15) in displacement coordinate systems at points GA and GB, respectively. 

(Real or blank) 


Rewirks : 1. Element identification numbers must be unique with respect to all other element 
identification nunbers. 

2. For an explanation of bar element geometry, see Section 1.3.2. 

3. Zero (0) must be used in fields 7 and 8 in order to override entries in these 
fields associated with F > ] in field 9 on a BAR0R card. 

4. If there are no pin flags or offsets, the continuation card may be omitted. 
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BULK DATA DECK 


Input Data Card CELAS1 Scalar Spring Connection 

Description : Defines a scalar spring element of the structural model. 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CELASl 

EID 

PIO 

G1 

Cl 

G2 

C2 



1 

CELASl 

2 

6 



8 

1 





J 


Field 


Contents 


EID 

PID 

G1, G2 
Cl, C2 


Unique element identification number (Integer >0) 

Identification number of a PELA:. property card (Default is EID) (Integer > 0) 
Geometric grid point identification number (Integer > 0) 

Component nunber (6 > Integer _ 0) 


R emarks .- 1. Scalar points may be used for G1 and/or G2 in which case the corre -.ponding Cl and'Or 
~ C2 must be zero or blank. Zero or blank may be used to indicate a groundeo'' 

terminal G1 or G2 with a corresponding blank or zero Cl or C2. If only scalar points 
and/or ground are involved, it is more efficient to use the CELAS3 card. 

2. Element identification nunbers must be unique with respect to aj_l_ other element 
identification nuitbers. 

3. The two connection points, (G1 , Cl) and (G2, C2), must be distinct. 


4. For a discussion of the scalar elements, see Section 5.6 of tne Theoretical 
fianual . 


* A grounded terminal is a scalar point or coordinate of a geometric grid point whose displacement 
is constrained to zero. 
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BULK DATA DECK 


Input Data Card C0NRHD Rod Element Property and Connection 

Description : Defines a rod element of the structural model without reference to a property cara. 
Format and Example : 


1 

2 

3 

4 


6 

7 

8 

9 

10 

C0NR0D 

EID 

Gl 

G2 

MID 

A 

J 

C 

NSM 


C0NR0D 

2 

16 

17 

23 

2.69 






Contents 

Unique element identification number (Integer > 0) 

Grid point identification nunbers of connection points (Integer ■> 0; -i r G2‘ 
Material identification number (Integer > 0) 

Area of rod (Real) 

Torsional constant (Real) 

Coefficient for torsional stress determination (Real) 

Nonstructural mass per unit length (Real) 


Field 

EID 

Gl, G2 

MID 

A 

J 

C 

NSM 


Remark s: 1. Element identification numbers must be unique with respect to other element 

” identification numbers. 

2. For structural problems, C0NR0O cards may only reference MATl material cards. 

3. For heat transfer problems, C0NR0D cards may only reference MAT4 or 
MATS material cards. 
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BULK DATA DECK 


Input Data Card CQUAD1 Quadrilateral Element Connection 

Description : Defines a quadrilateral membrane and bending element (QUADl) of the structural mode'. 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 ■ 9 

' c 

jcQUADl 

EID 

PID 

G1 

G2 

G3 

G4 1 

I • 

7H I 

1 

buAOl 

72 

13 

13 

14 

15 

iii 


1 


Field 


Con ten ts 


EID 

PID 

G1,G2,G3.G4 

TH 


Element identification number (integer >0) 

Identification number of a PQUADl property card (Default is EID) (Integer > 0. 

Grid point identification numbers of connection points (Inteoer > 0; 

G1 62 if G3 G4) 

Material property orientation angle in degrees iReal ) 

The sketch below gives the sign convention for TH. 



Remarks : 1. Element identification numbers must be unique with respect to aVl^ other element 

identification numbers. 

2. urid uoints G1 thru G4 must be ordered consecutively around the perimeter of the 
element. 

3. All interior angles must be less than 180°. 
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BULK DATA DECK 


Input Data Card CQUAD2 Quadrilateral Element Connection 

Description : Defines a homogeneous quadrilateral membrane and bending element (QUAD2) of the 
structural w)de1 . 



Field Contents 

EID Element Identification nunber (Integer > 0) 

PID Identification nutitoer of a PQUAD2 property card (Default is EID) (Integer > 0) 

G1,G2,G3,G4 Grid point Identification numbers of connection points (Integer > 0; 

G1 G2 ^ G3 ^ G4) 

TH Material property orientation angle in degrees (Real) ■ 

The sketch below gives the sign convention for TH. 



Remarks : 1. Element identification numbers must be unique with respect to ail, other element 

identification numbers. 

2. Grid points 61 thru G4 must be ordered consecutively around the perimeter of the 
element. 

3. All interior angles must be less than 180°. 
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BULK DATA DECK 


Input Data Card CSHEAR Shear Panel Element Connection 

Description : Defines a shear panel element (SHEAR) of the structural model. 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CSHEAR 

EID 

PID 

G1 

G2 

G3 

G4 




CSHEAR 

3 

6 

1 

5 

3 

7 





Field Contents 

EID Element identification number (Integer >0) 

PID Identification number of a PSHEAR property card (Default is EID) (Integer > 0) 

G1, G2, G3, G4 Grid point identification numbers of connection points (Integer > 0; 

61 F G2 G3 G4) 

Remarks : I. Element identification numbers must be unique with respect to alj_ other element 

identification numbers. 

2. Grid points G1 thru G4 must be ordered consecutively around the perimeter of the 
element. 

3. All interior angles must be less than 180". 
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BULK DATA DECK 


Input Data Card CTRIA1 Triangular Element Connection 

Description : Defines a triangular membrane and bending element (TRIAl) of the structural model. 



EID Element identification nunber (Integer > 0) 

PIO Identification nun*er of a PTRIAl property card (Default is EIO) (Integer > 0) 

G1,G2,G3 Grid point identification numbers of connection points (Integer > 0; 

G] f G2 f 63) 

TH Material property orientation angle in degrees (Real) - The sketch below gives 

the sign convention for TH. 



Remarks : 1. Element identification nunbers must be unique with respect to aVl_ other eleitent 

identification numbers. 

2. Interior angles must be less than 180°. 
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BULK DATA DECK 


Input Data Card CTRIA2 Triangular Element Connection 



Defines a triangular membrane and bending element (TRIAZ) of the structural model. 



Field Contents 

EID Element identification mjnber {Integer >0) 

PID Identification number of a PTRIA2 property card (Default is EID) (Integer > 0) 

G1,G2,G3 Grid point identification numbers of connection points (Inteoer > 0; 

G1 j< G2 ^ G3) 

TH Material property orientation angle in degrees (Real) - T^'e sketch below gives 

the sign convention for TH. 



Remarks : 1. Element identification numbers must be unique with respect to all_ other element 

identification nutibers. 

2. Interior angles must be less than 180°. 
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BULK DATA DECK 


Input Data Card F0RCE Static Load 

Description : Defines a static load at a grid point by specifying a vector. 
Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

in 

F0RCE 

SID 

G 

CID 

F 

N1 

N2 

N3 


i 

F0RCE 

2 

5 

6 

2.9 

0.0 

1.0 

0.0 




Contents 

Load set identification number (Integer >0) 

Grid point identification nunber (Integer > 0) 

Coordinate system identification number (Integer > 0) 

Scale factor (Real) 

Components of Vector neasured in coordinate svstein defined by CiD (Real; 
Nl^ + N2^ + N3^ > 0.0) 

Remarks: 1. The static load applied to grid point G is given by 

f i F fJ 

where N is the vector defined in fields 6, 7 and 8. 

2. Load sets must be selected in the Case Control Deck (L0AD=SID) to be 
used by NASTRAN. 


Field 

SID 

G 

CID 

F 

N1 ,N2,N3 


3. A CID of zero references the basic coordinate system. 
















BULK DATA DECK 


Input Data Card GRAY Gravity Vector 

Description : Used to define gravity vectors for use in determining gravity loading for the 

structural model . 



SID Set identification nunber (Integer > 0) 

CID Coordinate system identification number (Integer > 0) 

G Gravity vector scale factor (Real) 

Nl, N2, N3. Gravity vector components (Real; Nl^ + N2^ + N3^ > 0.0) 


Remarks : 1. The gravity vector is defined by 

g = G-(N1, N2. N3). 

2. A CID of zero references the basic coordinate system. 

3. Gravity loads may be combined with "simple loads" (e.g., F0RCE, M0MENT) only by 
specification on a L0AD card. That is, the SID on a GRAY card may not be the same 
as that on a simple load card. 

4. Load sets must be selected in the Case Control Deck (L0AD=SID) to be used by 
NASTRAN. 
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BULK DATA DECK 


Input Data Cara GRID Grid Point 


Description: Defines the location of a qeometric urid point of tne :.tructura1 model, the diroc- 

tTo'ns of fts displacement, and its permanent single-point constraints. 


> ormat and E xample : 


1 

2 

3 

4 

6 

6 

7 3 y 

Iv 

[grid 

ID 

CP 

XI 

X2 

i 

CO i .PS (>< 

j 

PID 

2 

3 

1.0 

2.0 

i 3.0 





Field 



)U,X2,X3 

CD 

PS 


Contents 

Grid point identification number (0<Integer<999999 ) 

Identification nunber of coordinate system in which the location of the grid 
point is defined (Integer > 0 or blank*). 

Location of the grid point in coordinate system C? (Real) 

Identification number of coordinate system in wh-ch displacements, degrees of freedom, 
constraints, and solution vectors are defined at the grid point (Integer > 0 or blank*) 

Permanent single-point constraints associated with grid point (any of the digits 
1-6 with no imbedded blanks) (Integer > 0 or blank*) 


R emarks ; 1. All grid point identification numbers must he unique with respect to all other 

structural, scalar, and fluid points. 

2. The meaning of XI , X2 and X3 depend on the type of coordinate system, CP, as 
follows: (see C0RO card descriptions) 


Type 

XI 

X2 

X3 

Rectangular 

X 

Y 

Z 

Cylindri cal 

R 

Q(degrees) 

Z 

Spherical 

R 

©(degrees) 

J(degrees) 


3. The collection of all CD coordinate systems defined on all GRID cards is called the 
Global Coordinate System. All degrees-of-freedom, constraints, and solution vectors 
are expressed in the Global Coordinate System. 


* See the GRDSET card for default options for fields 3, 7 and 8. 
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BULK DATA DECK 


Input Data Card MAT! Material Property Definition 

: Defines the material properties for linear, temperature-independent, isotropic 


Description 

materials. 



Field Contents 

MID Material identification nunber (Integer > 0) 

E Young's modulus (Real 2 0.0 or blank) 

G Shear modulus (Real >0.0 or blank) 

NU Poisson's ratio (-1 .0 < Real < 0.5 or blank) 

RH0 Mass density (Real) 

A Thermal expansion coefficient (Real) 

TREF Thermal expansion reference temperature (Real) 

GE Structural element damping coefficient (Real) 

ST, SC, SS Stress limits for tension, compression and shear (Real) (Required for Property 

Optimization calculations; otherwise optional if margins of safety are desired. V 


Remarks: 1. One of E or G must be positive (i.e., either E 0.0 or G > 0.0 or Doth £ and G .may 

he » 0.0). 

2 . If any one of E, G or NU is blank, it will be computed to satisfy the identity 
E = 2(1+NU)G; otherwise, values supplied by the user will be used. 

3. The material identification number must be unique for all MATl , MAT2 and MAT3 cards. 

4. MATl materials may be made temperature dependent by use of the MATT! card. 

5. The mass density, RH0, will be used to automatically compute mass for a’l structural 
elements except the two-dimensional bending only elements TRBSC, TRPLT and QDPLT. 

6. If E and NU or G and NU are both blank they will be both given the value 0.0. 

/. Weight density may be used in field 6 if the value - is entered on the ’’ARAM cars 
WTMASS, where g is the acceleration of gr.ivity. ^ 

8. Solid elements must not have NU equal to 0.5. 
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BULK DATA DfCK 



Field 

PID 

MIO 

A 

II, 12, 112 
J 

NSM 
Kl, K2 

Ci , Oi , Ei. Fi 


Contents 

Property identification number (Integer > 0) 
Material identification nunber (Integer > 0) 
Area of bar cross-section (Real) 

Area moments of inertia (Real, I^I^ 2 
Torsional constant (Real) 

Nonstructural mass per unit length (Real) 
Area factor for shear (Real ) 

Stress recovery coefficients (Real) 


Remarks : 1. For structural problems, P8AR cards may only reference MATl material cards. 

2. See Section 1.3.2 for a discussion of bar element geometry. 

3. For heat transfer problems, PBAR cards may only reference MAT4 or MATS 
material cards. 
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BULK DATA DECK 


Input Data Card PELAS Scalar Elastic Property 


Description : Used to define the stiffness, damping coefficient, and stress coefficient of a 
scalar elastic element (spring) by means of the CELASl or CELAS3 card. 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

3191 

PID 

K 

GE 

S 

PID 

K 

GE 

S 

1 

1 

PELAS 

7 

4.29 

0.06 

7.92 

27 

2.17 

0.0032 


r i 

1 



Contents 

Property identification nutiber (Integer > 0) 

Elastic property value (Real) 

Damping coefficient, g^ (Real) 

Stress coefficient (Real) 

Remarks : 1. The user is cautioned to be careful using negative spring values. (Values are 

defined directly on some of the CELASi card types.) 

2 . npo or two a'.astic spring properties may be defined on a single card. 

3. For a discussion of scalar elements, see 'lection 5.6 of the Theoretical Manual, 


Field 

PIO 

K 

GE 

S 
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BULK DATA DECK 


Input Data Card PQUAD1 General Quadri lateral Element Property 

Description : Defines the properties of a general quadrilateral element of the structural -"odel , 
including bending, membrane, and transverse shear effects. Referenced by the CQUADl card. 


Format and Example : 


1 2 3 4 5 6 7 8 ? 10 


Nuaoi 

PID 

MIDI 

T1 

MID2 

I 

MID3 

T3 

NSM 

abc 

POUADI 

32 

16 

2.98 


6.45 

16 

5.29 

6.32 

WXYZl 

Uc 

Zl 

Z2 








l+XYZl 

0.09 

-0.06 









Field 

PID 

MIDI 

T1 

MI02 

I 

MID3 
T3 
NSM 
Zl. Z2 

Remark' 


C ontents 

Property identification number (Integer 0) 

Material identification number for membrane (Integer 2 0) 

Membrane thickness (Real) 

Material identification number for bending (Integer 2 0} 

Area moment of inertia per unit width (Real) 

Material identification nunber for transverse shear (Integer ; 0) 

Transverse shear thickness (Real) 

Nonstructural mass per unit area (Real) 

Fiber distances for stress computation, positive according to the right-hand 
sequence defined on the CQUADl card (Real) 

1. All PQUADl cards must have unique property identification numbers. 

2. If T3 is zero, the element is assumed to be rigid in transverse shear. 

3. The membrane thickness, T1 , is used to compute the structural mass for this element. 
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BULK DATA DECK 


Input Data Card PQUAD2 Homogeneous Quadrilateral Property 

Decription : Defines the properties of a homogeneous quadrilateral element of the structural 
model, including bending, membrane and transverse shear effects. Referenced by the CQUAD2 card. 


Format and Example : 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

PQUAD2 

PID 

MID 

T 

NSM 

PID 

MID 

T 

NSM 



32 

16 

2.98 

9.0 

45 

16 

5.29 

6.32 



Field 

PID 

MID 

T 

NSM 


Contents 

Property identification nunber (Integer > 0) 
Material identification number (Integer > 0) 
Thickness (Real> 0.0) 

Nonstructual mass per unit area (Real) 


Remarks 


1. All PQUAD2 cards must have unique identification numbers. 

2. The thickness used to compute membrane and transverse shear properties is T. 

3. The area moment of inertia per unit width used to compute the bending stiffness is 
TV12. 

4. Cuter fiber distances of tT/2 are assumed. 

5. One or two homogeneous quadrilateral properties may be defined on a single card. 
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BULK DATA DECK 


Input Data Card PSHEAR Shear Panel Property 



Field Contents 

PIO Property identification number (Integer > 0) 

MID Material identification number (Integer > 0) 

T Thickness of shear panel (Real 0.0) 

NSM Nonstructural mass per unit area (Real) 

Remark s : 1. All PSHEAR cards must have unique identification numbers. 

2. PSHEAR cards may only reference MATl material cards. 

3. One or two shear panel properties may be defined on a single card. 
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BULK DATA DECK 


Input Data Card PTRIAI General Triangular Element Property 


Description : Defines the properties of a general triangular element of the structural model, 
including bending, membrane and transverse shear effects. Referenced by the CTRIAl card. 


Format and Example : 


1 2 3 4 5 ^ 7 8 9 10 


PTRIAI 

PID 

MIDI 

T1 

MID2 1 I 

M1D3 

T3 

NSM 

abc 1 

PTRIAI 

32 

16 

2.98 

9 j 6.45 

16 

5.29 

6.32 

QED ! 

♦be 


U 


r p . . .. 

i 

r 

i 

1 

— ,,J 


1 

1 


ED 


Contents 

Property identification number (Integer > C) 

Material identification number for membrane (Integer > 0) 

Membrane thickness (Real) 

Material identification number for bending (Integer >0) 

Area of moment of inertia per unit width (Real) 

Bending material identification number for transverse shear (Integer ; 0) 
Transverse shear thickness (Real) 

Nonstructural mass per unit area (Real) 

Fiber distances for stress calculations, oositive according to the right-hand 
sequence defined on the CTRIAl card (Real) 

Remarks : 1. All PTRIAI cards must have unique property identification numbers. 

2. If T3 is zero, the element is assumed to be rigid in transverse shear. 

3. The membrane thickness, T1 , is used to compute the structural mass for tnis element. 


Field 

PID 

MIDI 

T1 

MID2 

I 

MID3 

T3 

NSM 

/I. n 
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SULK DATA KCK 


Input Data Card PTRIA2 Hanogeneous Triangular Element Property 


Description : Defines the properties of a homogeneous triangular element of the structural model. 
Including membrane, bending and transverse shear effects. Referenced by the CTRIA2 card. 


Format and Example : 



Field Contents 

PID Property Identification number {Integer > 0) 

MID Material Identification number (Integer > 0) 

T Thickness (Real > 0.0) 

MSN Nonstructural mass per unit area (Real ) 

Hwarfes ; 1. All PTRIA2 cards must have unique Identification manbers. 

2. The thickness used to compute the membrane and transverse shear properties Is T. 

3. The area moment of Inertia per unit width used to compute the bending stiffness Is 
T*/12, 

4. Outer fiber distances of tT/2 are assumed- 

5. One or two homogeneous triangular element properties may be defined on a sinale 

card. ^ 
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BULK DATA DECK 


Input Data Card SPC1 Single-Point Constraint 

Description : Defines sets of single-point constraints. 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

a 

9 

i«j 

SPCl 

SID 

C 

Gl 

G2 

G3 

G4 

G5 

G6 

abc 

SPCl 

3 

2 

1 

3 

10 

9 

6 



ABC 


+bc 

G7 

G8 

G9 

-etc. - 



1 L. 

+BC 


8 

— 

J 




! i 


Alternate Form 


SPCl 

SID 

C 

GIDl 

"THRU" 1 GID2 



! 

i i 

SPCl j 

313 

12456 

6 

THRU 

32 1 


1 

1 

L, 3 


Contents 

Ident’ficaHen nu^’l.e'" of sinq’p-poif.t constrairit s=i. vinteger ■■ C 

Component numoer (Any unique combination of the digits ’-6 (with no i.~beaded 
blanks) when point identification numbers are grid ooints; must be null if oo- 
identification numbers are scalar points) 

Grid or scalar point identification numbers (Integer • 0) 


P^marks: T. Note tnat enforced displace.ments are not available via this card. As tranv continu.i 

tion cards as desired may appear when "THRU" is not used. 

2. A coordinate referenced on this card may net appear as a dependent coordinate '.r- a 

multipoint constraint relation, nor may it be referenced on a SPC, ji)Ml''l. 

SUPI3PT card. 

3. Sif'Uie-point constraint sets must be selected m the case Control Oerk i.S?C = clD' 
be used by NASTRAN. 

4. ‘.P( degrees of freedom may be redundantly specified as permanent constraints on V't 
GRID card. 

5. All grid points referenced by GIDl thru GID2 must exist. 


F ield 

SID 

C 

Gi , Gibi 
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BULK DATA DECK 


Input Data Card SPCADD Single-Point Constraint 


D escri pt ion : Defines a single-point constraint set as a union of single-point constraint sets 
defined via SPC or SPCl cards. 


Format and Example : 


1 2 3 4 5 6 7 8 9 10 


SPCADD 

SID 

SI 

S2 

S3 

S4 

S5 

S6 

S7 

abc 

SPCADD 

100 

3 

2 

9 

1 





+bc 

S8 

S9 

r 

-etc. - 


■ 

- — ~ 










J 




-etc. - 


F ield Conten ts 

Identification number for new sinnle-ooint constr.iint sot (Intcuer 0; f i01 or 
132 if sixisyttimetric) 

Si Identification numbers of single-poifit ^.onstrau.t bots via SfC cr SPCl 

cards (integer ■ 0; Slf. f Si 1 , 

A 

Fer.arks : 1. Single-point constraint sets must be selected in the Case Control ['eck (SPOSID) to 

be used by NASTRAN. 

2. No Si may be the identification number of a single-point constraint set defitied 
by at'other SPCADD card. 

3. The Si values must be unique. 

4. Set identification numbers of 101 or 102 cannot be used in a.v i symmetric problems. 
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Appendix C 
FAIRING STRUCTURE 




FAIRING NASTRAN MODEL 


Five-digit Con rod & bar element numbering scheme 


I tem 


Orientation 


In the X-direction 
In the Y-di rection 
in the Z-di rection 


10 

20 

30 


Grid No. 

XXX 

XXX 

XXX 


Four-digit CELASI element numbering scheme 


I tem 


Orientation 


Scalar spring in X-direction 1 
Scalar spring in Y-direction 2 
Scalar spring in Z-direction 3 


Grid No. 

XXX 

XXX 

XXX 


Tr iangu ] ar and quadrilateral e lament numbering scheme 

The smallest of the grid number used to define the element boundary is used 
as the identification number. 


Four-d i git t riangular Squadri lateral pro pert y identification numbering scheme 


I tem 


Matl 


Total thickness 


CQUADl, CTRIAl X 

CQUAD2, CTR1A2 X 


Y.YY 

.YYY 


Material numbering scheme 


Matl No. 

Matl 

1 

A1 umi num 

2 

Titanium 

3 

Steel 

4 

Fiberglass 

5 

Phenolic resin honeycomb 

6 

Aluminum honeycomb 
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Upper pivot fairings 




H23a 


H233 


-234 


-U235 


-J236 



Movab 1 e fairing 


Figure C-4. NASi'RAN ovenving fairing model - element identifications 























Figure C-6. - Bulkhead station 375.00 - pivot fairing. 


-,-*71 



local detail view from figure C-4. 




V i ew A Vi ev; 3 



View C V i ew 0 
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Figure C-8. - Undenving intermediate fairing panel - element 
identification, local panel views from figure C-3. 





Figure C-9. - Overwing fairing influence coefficient point location. 



201 209 



Figure C-10. - Underwing fairing influence coefficient point location. 
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r\BUi c;-i. - I'AiiuNC iM-rniiNci: i:oi;i-i-iL’u-:N'r points 


SIC 

point 

Description 

Coordinate 

X 

Y 

- 

H 

Upper pivot fairing 

924.600 

-150 . 500 

36.183 

B 

Upper pivot fairing 

1,015.250 

-130.500 

55.589 

B 

Upper pivot fairing 

895.200 

-146.500 

52.059 

14 

Upper pivot fairing 

924.600 

-146.500 

53.152 

17 

Upper pLvof fairing 

9(i9.n0 

-145.000 

55.2"5 

20 

Upper pivot fairing 

1,01.5.250 

-14~.000 

55.100 

24 

Upper pivot fairing 

924.600 

-165.00 

29.11" 

30 

Upper pivot fairing 

1,013.250 

-165.000 

31.48" 

44 

Lower pivot fairing 

910.500 

-153.500 

-5.068 

49 

Lower pivot fairing 

910.500 

-147.291 

-5.282 

54 

Lower pivot fairing 

970.800 

-163.975 

-3.189 

57 

Lower pivot fairing 

1,012.100 

-174.784 

-1.594 

58 

Lower pivot fairing 

883.800 

-154.600 

-1.412 

60 

Lower pivot fairing 

910.700 

-166.237 

-1.307 

62 

Loiver pivot fairing 

957 . 500 

-169.653 

-2.092 

64 

Lower pivot fairing 

959.250 

-171.982 

-5.650 

65 

Lower pivot fairing 

9^0.800 

-173.116 

-5. Sob 

66 

Lower pivot fairing 

983.250 

-174. 538 

-2.756 

70 

Lower pivot fairing 

912.250 

-181.000 

1.196 

75 

Lower pivot fairing 

970.800 

-182.158 

-1.238 

78 

Lower pivot fairing 

1,012.100 

-185.769 

-0.913 

85 

Overwing movable fairing 

1,133.579 

-142.022 

16.105 

96 

Overwing movable fairing 

1,111.821 

-152.848 

25.367 

112 

Overwing movable fairing 

1,069.687 

-155. ~64 

50.181 
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TABLE C-I. - FAIRING INFLUENCE COEFFICIENT POINTS - Continued 


SIC 

point 

Description 

Coordinate 

X 

Y 

- 

114 

Ovenving Movable Fairing 

1,077.685 

-169.854 

27.968 

116 

Overwing inovalile fairing 

1 ,084.236 

-185.004 

24.879 

127 

Overwing movable fairing 

1,04~.467 

-184.869 

27.621 

128 

Forward intermediate fixed fairing 

1,036.000 

-165.480 

51 . 155 

133 

Forward intermediate fixed fairing 

1,043.500 

-124.250 

52.240 

134 

Forward intermediate fixed fairing 

1,065 . 500 

-142.580 

30.52" 

135 

Forward intermediate fixed fairing 

1,065.500 

-124.230 

50.415 

136 

Forward intermediate fixed fairing 

1,065.500 

-110.360 

30.988 

137 

Forward intermediate fixed fairing 

1,085.100 

-124.230 

27.691 

140 

Forward intermediate fixed fairing 

1,096.100 

-115.000 

28 . 336 

145 

Upper pivot fairing 

924.600 

-187.500 

19.5"" 

151 

Upper pivot fairing 

968.300 

-18". 500 

24.006 

157 

Upper pivot fairing 

1,011.000 

-18". 500 

24.580 

201 

Underwing fairing interned i ate panel 

1 , 165 . 87 

-116.50 

4.9520 

203 

Underwing fairing intermediate panel 

1,165.87 

-152.20 

4.5284 

205 

Underwing fairing intermediate panel 

1,165.87 

-151. “1 

3 . 5556 

207 

Underwing fairing intermediate panel 

1,165.87 

-160.86 

3.1902 

209 

Underwing fairing aft panel 

1,214.50 

-117.10 

5.6481 

211 

Underwing fairing aft panel 

1,215.10 

-135.91 

4.9154 

213 

Underwing fairing aft panel 

1,212.29 

-143.76 

4.4829 

215 

Underwing fairing aft panel 

1,211.16 

-157.56 

5 . 8886 

217 

Underwing fairing intermediate panel 

1,129.20 

-175.46 

1.8669 

219 

Underv\fing fairing intermediate panel 

1,147.54 

-173.61 

n 

221 

Undenting fairing intermediate panel 

1,165.87 

-171. "6 

H 
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TABLE C-I. - FAIRING INFLUENCE COEFFICIENT POINTS - Concluded 


SIC 

point 

Description 

Coordinate 

X 

Y 

1 1 

223 

Underwing fairing intermediate panel 

1,180.35 

-170.50 

5.1086 

225 

Underwing fairing intermediate panel 

1,190.90 

-169.24 

5.5646 

111 

Underwing fairing aft panel 

1,191.00 

-169.24 

5.5646 

229 

Underwing fairing aft panel 

1,210.56 

-166.97 

5.4685 

231 

Undenting fairing aft panel 

1,227.19 

-165.00 

5.5546 

233 

Underwing fairing aft panel 

1,245.50 

-165.10 

5.6562 
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Vertical deflection 



130 UO 150 160 .170 180 

Fuseiage Y-coordinates - in. 

Figure C-11. - Deflections for XASTRAN upper pivot fairing for unit 
loads applied at SIC 4 and 14. 







Vertical deflections 



Figure C-13. - Deflections for NASTRAN undenving intermediate fairing 
for unit loads applied at SIC 205 and 205. 


69 



Vertical deflections 





Vertical deflection 






Airloads Research Study - Fairing Stmcture 





1 h 

r BULK 

r A . A 5 r H 0 



CUt.j 








1 M 

• 1 • • 2 

.. 3 

.. 5 

• • 5 • • 

8 .. 7 .. 6 

• • 

9 .. 10 . 

2- 

C0A' IfilOl 

19992 

lu 1 

111 

10.1 

1 

baric 101 

3- 

♦ AtUulCl 56 

1 







CFAv 10111 

19992 

111 

121 

10*0 

1 

E7P10111 

B- 

♦AkIOIII 

1 






6- 

CRA7 1J119 

11531 

llA 

124 

lO.J 

1 


T- 

C6A< 10121 

19992 

121 

131 

10.0 

1 

BAR10121 

8- 

♦ A.nU21 
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n- 

CbH lul29 

13 c 30 

12A 

134 

10.0 

1 


10- 

C6Ar 10131 

19992 

1^1 

141 

IP.O 

1 

EA.M0131 

11- 

♦ A;.,l')131 56 

1 






12- 

CbA^ 10139 

19 66 1 

13«- 

144 

1 J.O 

1 


11- 

CPA? 10191 

1 9992 

111 

151 

10.0 

1 

EAR 10 141 

14- 

♦ AKltll91 

1 56 






15- 

C8A7 10199 

1993'J 

191 

154 

10. 0 

1 


16- 

C6A« 1C159 

16100 

15t 

164 

10.1 

1 


1/- 

CBA< 1C201 

19992 

cC 1 

211 

10. 0 

1 

EAR10201 
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